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Abstract 

The time evolution of a wormhole in a Friedmann universe approaching the Big 
Rip is studied. The wormhole is modeled by a thin spherical shell accreting the 
superquintessence fluid - two different models are presented. Contrary to recent 
claims that the wormhole overtakes the expansion of the universe and engulfs it 
before the Big Rip is reached, it is found that the wormhole becomes asymptotically 
comoving with the cosmic fluid and the future evolution of the universe is fully 
causal. 
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1 Introduction 


The 1998 discovery that the expansion of the universe is accelerated |I|, confirmed by the 
analysis of the cosmic microwave background power spectrum im has led to postulate, 
as a possible explanation within the context of general relativity, the existence of a form 
of dark energy or quintessence with negative pressure P < —p/3 (where p is the dark 
energy density). The dark energy is dynamically irrelevant in the earlier stages of the 
evolution of the universe and starts dominating the cosmic dynamics only at recent times. 
At present the observational data favour an even more exotic dark energy with effective 
equation of state P < —p violating the weak energy condition 0) . This violation is rather 
disturbing because it allows, in principle, for exotic solutions of general relativity such 
as wormholes and warp drives and the possibility of time travel associated with them. 
Another concern is that a universe dominated by dark energy with effective equation of 
state parameter w = P/p < — 1 may end its existence at a Big Rip singularity at which 
the scale factor and the energy density and pressure of quintessence diverge at a finite 
time in the future |5]. Dark energy with w < — 1 is called superquintessence or phantom 
energy. 

Recently, the evolution of a wormhole embedded in a Friedmann-Lemaitre-Robertson- 
Walker (hereafter “FLRW”) universe approaching the Big Rip was studied by Gonzalez- 
Diaz [ 0 ], with the conclusion that the wormhole accreting superquintessence expands 
faster than the background FLRW universe and that the radius of the wormhole throat 
diverges before the Big Rip is reached - thus the wormhole engulfs the entire universe, 
which will re-appear from the other wormhole throat. The resulting spacetime is not 
globally hyperbolic and is acausal with closed timelike curves threading the wormhole 
throat. Such bizarre scenarios are potentially of interest as constraints: if it can be 
established that phantom energy leads in principle to unacceptable consequences, this 
may be sufficient to rule out its existence. 

The analysis of Ref. is based on a qualitative estimate of the accretion rate of 
phantom energy onto the wormhole and the rate of variation of the throat radius, which 
borrows from a recent study of a similar problem of accretion of dark energy onto 
black holes [?|. The purpose of this work is to present exact solution models of a 
wormhole immersed in a spatially flat FLRW universe and to compare the expansion of 
the wormhole throat with that of the cosmic fluid as the Big Rip is approached. Two 
different and rather general wormhole models are studied. In the first model it turns 
out that, in the simplifying approximation of stationary accretion advocated in Ref. (ij. 
the wormhole is asymptotically comoving with the FLRW background as the Big Rip is 
approached - the size of the wormhole with respect to a comoving ruler does not increase 
and the universe cannot disappear within the wormhole. This conclusion holds also in 
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the second wormhole model (presented in Sec. 5), which is more general. 


2 A thin shell wormhole embedded in a spatially flat 
FLRW universe 


We study an exact solution of the Einstein equations describing a spherically symmetric 
thin shell wormhole embedded in a spatially flat FLRW universe described by the metric 




ds 2 = —dt 2 + a 2 (t ) 


dr 2 + r 2 


(d6 2 


+ sin 2 6 d </k 


( 2 . 1 ) 


in comoving polar coordinates ( t , r, 6, p). The thin shell is located on a surface E of con¬ 
stant comoving radius. Two FLRW regions joining smoothly on the surface E constitute 
the regions “above” and “below” E. The assumption that the line element is given by 
eq. EH) implies that the wormhole shell does not perturb the surrounding universe (see 
the discussion in Sec. 4). The equation of the shell E is 


r = 


a(t) 


( 2 . 2 ) 


where R(t) is the comoving radius of the shell and the function a{t) is introduced for 
later convenience. The normal vector to the shell, obtained by differentiating eq. EI21 

is 

N a = (-a t e- Q ,cr 2 ,0,0) (2.3) 

where a subscript t denotes differentiation with respect to the comoving time t of the 
FLRW background. The norm squared of N a is 

g a0 N“N^ = -af e~ 2 “ + \=(£) , (2.4) 

a z \a J 

where 

/3 = \J 1 — a 2 a 2 e~ 2ot = y/l — v 2 (2.5) 

can be interpreted as the inverse of the Lorentz factor constructed with the velocity of 
the shell relative to the background 


v = —a t R . 


By normalizing N a one obtains the unit normal to the surface E 


n 


r — 


(~a t R 

\ 13 



( 2 . 6 ) 


(2.7) 
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and the three-metric on the surface £ is given by 


ds 2 


s 


—/3 2 dt 2 + R 2 (t) (d6 2 + sin 2 0 dip 2 ') 


—dr 2 + R 2 (r)dQ 2 


( 2 . 8 ) 


in coordinates (t, 9, ip) on £, where r is the proper time in the shell’s frame (defined by 
dr = j3dt ) and dQ 2 is the line element on the unit two-sphere. By using the tetrad 

b? t >’ e TO’ e w} = {M“A“ K} (2-9) 


(a, b = t,9,<p), one computes the extrinsic curvature of £ 

K aP = e<“> e f W a n b = e £> ef (, d a n b - T> c ) . 
Its components are found to be 


KK 


~r 3 ‘ 


'a t R s 




and its trace is 


K ° 0 


K\ 


(5 at Rt 

R /T~ 


K 


K t t + K e e + K\ = -^d t 



2/3 2ott Rt 

+ ~R — 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 


Since the shell has no interior and two FLRW regions match on E, the jumps of A" 1 ’,, 
and of K at E, which appear in the Einstein equations at the shell, are given by 



[. K\} = 2 K a b , [K] = 2 K 

(2.15) 

The Einstein equations at £ are ’2] 



[K\ - K5\\ = -8irS a b , 

(2.16) 

where 

Sab = (a + P( E)) ] + P (s) g ab | s 

(2.17) 

is the stress-energy tensor of the (exotic) matter living on the shell and u ^ is the four- 
velocity of the shell, cr and P(s) are, respectively, the (surface) energy density and the 
pressure on the shell. The time-time component of the Einstein equations yields 


1 f P a t R t \ 

2tt l R p ) ’ 

(2.18) 
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while the 99 or the (pip component yields 


P, 


(£) 


1 [ 1 rj (ot t R\ (3 a t Rt 

An [ (3 2 1 ^ (3 J + R (3 


a _ d t ( u a n a ) 
2 An(3 2 


(2.19) 


where u M is the four-velocity field of the cosmic fluid. The effective equation of state of 
the exotic matter on the shell is given by 


<7 + 2 P(s) — 


1 . d t ' atR 


2n(3 2 V P J 

while the material energy residing on the shell is 

M = AnR 2 a = -2 ((3R 


2n(3‘‘ 


a t R t R 2 


dt (u a n r 


(3 


( 2 . 20 ) 


( 2 . 21 ) 


The relative speed v = —a t ae~ a between the shell and the background fluid, [3 = 
a /1 — v 2 , and u a n a = —v//3, are constant when the matter on the shell satisfies the 
equation of state P^) = —cr/2. 

The Einstein equation 


where 


k ab S ab =[T aP n a n 0 } , 

(2.22) 

zy(+) , jy(-) 

K ab = ab + ab =0 
ab 2 

(2.23) 


is the mean of the extrinsic curvature on both sides of the shell, provides the rather 
obvious matching condition for the energy density and pressure of the FLRW cosmic 
fluid on both sides of E: 


p (+) = p (-) f P (+) = p(-) . (2.24) 

ft is sometimes convenient to use the parameter y defined by 

v = —a t R = —d[3R = tanh y , (2.25) 

or by coshy = /3 -1 , where an overdot denotes differentiation with respect to the comov¬ 
ing time of the shell r and f — ft/P for any differentiable function /. 

The four-velocity of the shell 

r\ r r 0i 1 dr a 1 / x 

(2 - 26) 
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has the properties 


U (H) n ot — 0 , 


W(S) M « S) = _1 ' 


(2.27) 

(2.28) 

(2.29) 


The conservation equation (eq. (5) of Ref. [5j) projected along the four-velocity of the 
shell yields 


SaVfo = - 


U 


(S) 




(2.30) 


The left hand side S a b . b u^ = ~S b & of eq. (12.3011 reduces to — (au^^j — 

P(S)W(s) b i where (= M /A, A = AnR 2 (t) is the (proper) surface area of the 
shell, and M = cr.A is the material energy located on the shell. Similarly, it is found 
that u l {f = A/A. The stress-energy tensor of the cosmic perfect fluid in the FLRW 
background is 

T af) — (P + p) u a up + Pg a p , ( 2 - 31 ) 

where u a is the four-velocity of the FLRW comoving observers. Then 

T$n 0 = (P + p) (u^u a ) (yPnp) + Pu^ n Q . ( 2 . 32 ) 

By using the relations u%^u a = — /3 _1 , u^np = aR, and n" E yn Q = 0 one obtains 

( 2 . 33 ) 


u 


^np=-^j^-aR. 


Since the unit normal ip to £ has opposite sign “above” and “below” £, the jump in 
this quantity is 


u i s) TSn* 1 = - - (P + p) aR 


A- i i3 

and the conservation equation yields, using aR = —v/(3 , 

2 

W 2 


M + P{t.)A = —(P + p)Av 


(2.34) 


(2.35) 


This equation describes the rate of accretion of the cosmic fluid by the wormhole and 
can be interpreted as follows. The quantity (P + p) v on the right hand side is the flux 
density of the cosmic fluid crossing the shell £ radially, the factor 2 arises because the 
outflow is from both faces of the shell, one factor /P 1 conies from the relativistic mass 
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dilation, while another factor f3 comes from Lorentz contraction in the radial direction 
due to the relative motion of the shell and the FLRW background. Note that in a de 
Sitter background enjoying the equation of state P = —p there is no accretion on the 
shell and static solutions with both M and R vanishing (such as those considered in 
Ref. |TDJ or their generalizations) become possible. 

If the strong energy condition is satisfied by the cosmic fluid and P + p > 0 then a 
wormhole shell expanding relative to the cosmic substratum (v > 0) accretes positive 
cosmic fluid energy while a shell contracting (v < 0) relative to the cosmic substratum 
experiences an energy outflow through £. 

3 v =constant solutions 

In the special case in which the relative radial velocity of the cosmic fluid and the 
wormhole shell is constant, 


v = — a t R = tanh y = u 0 , (3.1) 

the equation of state of the exotic matter on the shell is P(^) = —cr/2 and one can 
eliminate R(t) betwen eqs. (13.11) and m obtaining 

- a(t) d t (e" Q(t) ) + v 0 = 0 . (3.2) 

We are interested in a FLRW universe approaching the Big Rip and therefore, as done 
in Ref. |5], we assume for simplicity that the equation of state of the cosmic fluid is 
constant, P = wp with w < —1. This guarantees the occurrence of a Big Rip described 
by the form of the scale factor 


a(t) = a 0 ( t r i P ~ t) 3( “ +1) , (3.3) 

which diverges together with the energy density p = po a 3< - u,+1 ' ) and the pressure P = wp 
as t — > t r i p . Under this assumption the constant v solution of eq. (EH for the motion of 
the wormhole shell is given by 

e -“W = Cv o (trip - t )rar + e-“° , (3.4) 

where C and «o are constants. The comoving radius of the shell is 

R(t) = a(t ) e“ a(t) ~ a(t) e~ ao : (3.5) 
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it scales asymptotically like the scale factor a(t) and hence the wormhole does not 
overtake the expansion of the universe, contrary to what is suggested in Ref. [0]. 

The fact that the shell cannot expand faster than the universe indefinitely can also 
be seen by differentiating eq. Q) with respect to t, which yields 


Rt_ 

R 


v 

h + r' 


(3.6) 


where H = a t /a is the Hubble parameter of the FLRW universe. When v is constant and 
the radius of the shell expands to infinity the term v/R in eq. (18.611 becomes negligible 
and the expansion rate of the shell (with respect to comoving time) coincides with the 
expansion rate of the background universe, even if the shell starts out expanding at a 
faster rate than the background. The final state is indistinguishable from one with v = 0 
and it is the same irrespective of the initial conditions. (Strictly speaking, this argument 
does not assume that the Big Rip is approached, but only that the universe expands 
forever and v =const., which in turn implies that the wormhole shell expands to infinity, 
R(t) —> Too). 

Perhaps a more physical way of looking at this aspect is the following: if the shell 
were to expand faster than the cosmological background in the approach to the Big Rip, 
and given that H —> Too near the Big Rip, the relative speed v between the shell and 
the cosmic fluid would increase its magnitude without bound. This would contradict 
the fact that |u| must be bound by unity, or else the wormhole shell becomes tachyonic. 
Hence the asymptote |u| ~ 1 must be approached before the Big Rip is reached, with 
dtv —> 0. A comparison of these results with those of Ref. |Bj is given in Sec. 6. 


4 Gravitating mass of the wormhole 

We briefly discuss our assumption that the wormhole considered here does not perturb 
the FLRW surroundings. This assumption must correspond to a zero gravitational mass 
for the wormhole. In the static limit a = 1 the mass-energy on the shell is the only 
mass-energy in the entire spacetime. Because the latter is asymptotically flat, the total 
gravitational mass is given by the Tolman mass EU; for a shell with surface stress- 
energy tensor S a b in an asymptotically flat spacetime the Tolman mass is given by the 
expression EO 

M t = JjE^m(-S° 0 +S 2 2 +Sl) 

= fj dip j* ddR 2 (3sinO(a + 2P^) = ~ • ( 4 - 1 ) 
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If a t R is constant, corresponding to constant v and /?, the Tolrnan mass Mt vanishes 
identically and the wormhole does not affect the Minkowski background both above 
and below £. Note that not only the material mass but also the pressure P(s) 

contributes to the Tolrnan mass, and that the pressure is adjusted to compensate the 
contribution of the material mass in such a way that the Minkowski background is not 
altered. 

For a wormhole embedded in a FLRW space the Tolrnan mass is replaced by the 
Hawking quasi-local mass mm- since the metric m outside the wormhole shell is ex¬ 
actly a FLRW one, the quasi-local mass reduces to that of a FLRW space, E = 47rr 3 p/3, 
with no contribution from the shell. One concludes that the gravitational mass of the 
wormhole is zero. If instead a spherical black hole is embedded in a FLRW background, 
as described e.g., by the McVittie solution na, the quasi-local mass coincides with the 
black hole mass times the scale factor a(t) fH]. As an astrophysical object, a zero-mass 
wormhole construct located in a galaxy would not perturb the orbits of stars that do 
not fall directly into its throat, and it would not cause gravitational lcnsing (lensing by 
negative mass wormholes is studied in Refs. dl). 


5 More general wormholes 


The wormhole model considered in sections 3 and 4 suffers from the intrinsic limitation 
that the wormhole shell is adjusted so that it does not perturb the surrounding cosmo¬ 
logical background. It is interesting to ask whether removing this assumption has an 
effect on the results of Sec. 3 about the relative motion of the shell and the background 
fluid as the Big Rip singularity is approached. To this purpose, consider the McVittie 
metric m 


ds 2 


B(t,r ) 
A(t,r) 


dt 2 + a 2 {t)A A (t, r) (dr 2 + r 2 dfl 2 ^j 


(5.1) 


in isotropic coordinates (t, r, 6, ip), where 


A(t,r) = 1 + 



B(t,r ) 


771 (t) 

2 r 


(5.2) 


The modification of the FLRW metric introduced by a non-vanishing function m(t) is 
caused by a spherical wormhole shell £ located at the radius 

m 

a(t)A 2 (f,r s ) ' 
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r = r^it) 


(5.3) 








The three-dimensional metric on E is given by 


ds 2 L = 


B 

A 


sechy) dt 2 + R 2 {t ) dQ 2 


in coordinates (t, 9, p) on E and with the function y defined by 

A 3 (f,r s ) dr s 


tanhy(i) = a(t) 


B (t, r s ) dt 


The unit normal to E is 


n u = A 1,0,0 


dt 


X = m4 2 cosh y , 


the four-velocity of the shell is 


W (S) = ^ cosh X ^1, 0, o) , 


and the four-velocity of the cosmic fluid is 

(A 


u»= (|,0,0,o) 


The projection of the shell four-velocity on the cosmic fluid four-velocity is 

1 


“(s) “t> = — cosh x = 


VT - 


(5.4) 


(5.5) 


(5.6) 


(5.7) 


(5.8) 


(5.9) 


where the three-dimensional velocity of the shell relative to the cosmic fluid is defined 


as 


v = tanhy(f) . (5.10) 

By differentiating ( R/a ) and using the expressions (|5.2(1 of A(t, r ) and B(t , r) one obtains 
an equation regulating the dynamics of the wormhole shell: 


dry 


dt AB 


d f R\ ^ 171 

dt \ a J dt 


(5.11) 


The extrinsic curvature of the shell E given by eq. S has the only nonzero compo¬ 
nents 

K r. _ \ 1 dr s dm\ m 1 

* B 2 ^ dt 2 ABrn dt dt) a 2 A 5 Brl\ ’ 1 J 
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I<°n = = A 


B 


+ 


A 2 dry 


H + 


1 dm 


(5.13) 


a 2 A 5 ry B 2 dt ' ir s dt 

where, again, H = a t /a. The Einstein equations ( 12 .Kill at £ then yield the surface 
density a and the pressure P(£) of the shell material 


K e 


a = 


e 


2n 

a + 2 P(s) = 


1 

2n 


B 


aA 3 ry 


/ 

cosh x + — \ H + 


B 


1 dm 
Ar E dt 


sinh x 


K\ 

2vr 


A 

27T 


A 2 (d 2 ry 1 


dry dm N 


B 2 \ dt 2 ABry dt dt 


+ 


m 


a 2 r^A 5 B 


(5.14) 


(5.15) 


and the mass of exotic matter on the wormhole shell can be written, using eq. (15.31) . 


M = 4vr R 2 {t)a = -2R 


B , AR 
-co S h x + — 


H + 


1 dm 


Ary dt 

The cosmic imperfect fluid is described by the stress-energy tensor 

T a p (P T P) ^a'U'P T PPa/3 T Qa^P T QpU a , 


sinh x 


(5.16) 


(5.17) 


where the purely spatial vector (with q^u^, = 0) is the radial energy flux density. The 
covariant conservation equation of mass-energy for the shell, eq. (ism yields 


1 dM 


A \dry 


+ P (£) 


dA 

dry 


= 2A cosh y < (P + p) 


A dry 
B dt 


A 2 (dry 


B 2 \ dt 


a 2 A* + 1 


(5.18) 


where A = 47 rP 2 (f) and Ty is the proper time of the shell defined by the three-metric 
(15.41) as 

dry — — seedy dt . (5.19) 

The Einstein equation Gl = 87tT ( ) = 87 \q r u t determines the radial flux density q r as 


2m 
AB rf 


H + 


1 dm\ 


= 8tt q r ut , 


m dt J 

where rs(t,r ) = a(t)A 2 (t,r)r. Further use of eqs. (15.61) and m yields 


(5.20) 


dM dA 2 
dt ^ dt (3 


B A (1 + v 2 ) d ( ma ) 

- {P + P )A V + - ^ I 


(5.21) 
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or 


dM dA 

+ %> 


dr 


dr y 


2 

0 




d ( ma ) 

dry. 


(5.22) 


in terms of te, where /3 = \/l — v 2 again. Eq. 1 IT 22 I 1 reduces to eq. 'ITTli when m{t) = 0. 
The second term on the right hand side of eq. (EQ) is present even when v = 0 and 
describes a contribution to accretion onto the shell due to the radial energy flux of cosmic 
fluid - this term vanishes if 

~ = ~H, (5.23) 

m dt 

equivalent to m{t) = const. /a(t), or q r = 0 and Gj = 0. The quantity m(t ) a{t ) appear¬ 
ing in the second term on the right hand side of eq. (EQ) coincides with the Hawking 
quasi-local mass j2T]. The “Schwarzschild mass function” ms ( t,r ) instead is defined by 


1 - 


2m s (t, r, _ p 


= 9 ap (V a r s ) (Vpr s ) . 


(5.24) 


r s (t, r) 

The dehnitions of A, B, and rs and the time-time component of the Einstein equations 


yield 


A\ 2 ( TT 1 dm\ 2 

3I b) [ h + 7a7f) = 8 ^ 


ms ( t , r) = m(t) a(t ) +r 3 s p 


(5.25) 


(5.26) 


which reduces to the usual Schwarzschild mass in the absence of cosmic fluid. 

We are left with the problem of solving for the dynamics of the wormhole shell, which 
is determined by eq. (I5TTT1) . The latter can be written as 


d / R\ a d‘ m v f 2 

dt \ a J dt a V A, 


(5.27) 


During the approach to a Big Rip singularity in which the scale factor a(t) has the form 
(D, the metric component goo = — {B/A) 2 stays finite even if m(t) diverges. It seems 
reasonable to require that both A and B be finite and that the divergence be contained 
in a(t), keeping in mind that the mass of the wormhole shell (which is allowed to diverge) 
is not m(t), but is obtained by integrating eq. (15. 2 lib On the other hand, if v is finite 
(v < 1) eq. (15.271) reduces to the asymptotic equation 


d (R\ . dm 
— — )-A— = 0 
dt \ a J dt 


(5.28) 
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By using the implicit definition (15.31) of ry one obtains B ( t,rdr^/dt = 0, which 
has the constant solution = C , or R(t) = a(t) C + m(t) + rn 4 ^ . Again, approaching 
the Big Rip, the wormhole shell becomes comoving with the cosmic substratum. 


6 Discussion and conclusions 


We are now ready to discuss the difference between our results and those of Gonzalez- 
Diaz p. Gonzalez-Diaz begins by considering a static wormhole solution appearing in 
the appendix of Ref. [ IS] (eqs. (A.28) of [IHJ) and consisting of a spherical thin shell of 
exotic matter with throat radius b 0 and mass /i related by .TS] 


R = 


nbo 

2 


( 6 . 1 ) 


Gonzalez-Diaz extrapolates eq. eb to a time-dependent wormhole embedded in a 
FLRW universe and proceeds to estimate the accretion rate on such a wormhole by 
adopting the formula derived for stationary accretion onto a black hole embedded in a 
FLRW universe jl] 

/i = AnDp 2 (P + p) , (6.2) 

where D is a constant. By combining eqs. (EH) and E2P, Gonzalez-Diaz obtains an 
equation for the rate of change of the wormhole throat b 0 and proceeds to solve it 
in a FLRW universe approaching the Big Rip, concluding that the wormhole shell of 
radius b 0 (t ) expands faster than the universe and ends up engulfing it and destroying 
global hyperbolicity. Although treating the wormhole as a black hole with negative mass 
seems reasonable, unfortunately eq. EH) does not hold for a time-dependent wormhole 
embedded in a FLRW universe - cf. eq. (ESP or eq. (IT2TD for the first wormhole 
model, or eq. (15.161) for the second wormhole model presented. The time evolution of 
the wormhole shell cannot be guessed a priori but needs to be derived from detailed 
models like the ones presented here. 

The first wormhole model studied in this paper is appropriate for elucidating the 
physical situation considered in Ref. [S]. In fact, by comparing eq. (2) of Ref. pa] and our 
eq. (EQ51) . one sees that the regime of stationary accretion onto the wormhole considered 
in Ref. L 6] corresponds to the situation v =constant in our formalism, for which the 
dynamics of the shell has been explicitly solved (the static limit of such a solution would 
correspond to a = 1 and M = SttP^RR). In this case the wormhole shell ends up 
expanding at the same rate as the universe in which it is embedded. However, this 
wormhole model has one limitation, namely the assumption that the wormhole does 
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not perturb the surrounding FLRW universe. More general exact solutions describing 
a wormhole embedded in, and modifying, the surrounding cosmological background 
are presented in Sec. 5. This more general class of exact solutions contains three free 
functions a(t), m(t), and R(t). By imposing the form (Id.hi) of the scale factor a(t ) 
appropriate to the description of a Big Rip and by noting that the metric component 
g 00 = — ( B/A ) 2 is bounded even if m(t) diverges, one deduces again that the wormhole 
shell becomes comoving with the cosmological background as the Big Rip is approached. 


Acknowledgments 

This work was supported by the Natural Sciences and Engineering Research Council of 
Canada ( NSERC ). 


References 

[f] A.G. Riess et al, Astron. J. 116, 1009 (1998); S. Perlmutter et ai, Nature 391, 51 (1998); 
A.G. Riess et a/., Astron. J. 118, 2668 (1999); S. Perlmutter et ai, Astrophys. J. 517, 
565 (1999); A.G. Riess et al, Astrophys. J. 536, 62 (2000); A.G. Riess et al ., Astrophys. 
J. 560, 49 (2001). 

[2] C.B. Netterfield et al ., Astrophys. J. 571, 604 (2002). 

[3] G. Hinshaw et al, Astrophys. J. Suppl. 148, 135 (2003). 

[4] A. Melchiorri, L. Mersini, C.J. Odrnan and M. Trodden, Phys. Rev. D 68, 043509 (2003); 

R. R. Caldwell, Phys. Lett. B 545, 23 (2002). 

[5] R.R. Caldwell, M. Kamionkowski and N.N. Weinberg, Phys. Rev. Lett. 91, 071301 (2003); 

S. M. Carroll, M. Hoffman and M. Trodden, Phys. Rev. D 68, 023509 (2004). 

[6] P.F. Gonzalez-Diaz, Phys. Rev. Lett. 93, 071301 (2004). 

[7] E. Babichev, V. Dokuchaev and Yu. Eroshenko, Phys. Rev. Lett. 93, 021102 (2004). 

[8] We use units in which the speed of light c and Newton’s constant G assume the value 

unity, the metric signature is — and the other notations follow those of Ref. fSjj- 

[9] C. Barrabes and W. Israel, Phys. Rev. D 43, 1129 (1991). 

[10] J.P.S. Lemos, F.S.N. Lobo and S. Quinet de Oliveira, Phys. Rev. D 68, 064004 (2003). 


13 


[11] R.C. Tolman, Relativity, Thermodynamics, and Cosmology (Clarendon Press, Oxford, 
1962). 

[12] V. de La Cruz and W. Israel, Phys. Rev. 170, 1187 (1968). 

[13] S.W. Hawking, J. Math. Phys. 9, 598 (1968). 

[14] S.A. Hayward, Phys. Rev. D 49, 831 (1994). 

[15] G.C. McVittie, Mon. Not. R. Astr. Soc. (Lon.) 93, 325 (1933). 

[16] B.C. Nolan, Phys. Rev. D 58, 064006 (1998). 

[17] J.G. Cramer, R.L. Forward, M.S. Morris, M. Visser, G. Benford and G.A. Landis, Phys. 
Rev. D 51, 3117 (1995); D.F. Torres, G.E. Romero and L.A. Anchordoqui, Phys. Rev. D 
58, 123001 (1998); Mod. Phys. Lett. A 13, 1575 (1998); L.A. Anchordoqui, G.E. Romero, 
D.F. Torres and I. Andruchow, Mod. Phys. Lett. A 14, 791 (1999); M. Safonova, D.F. 
Torres and G.E. Romero, Mod. Phys. Lett. A 16, 153 (2001); E. Eiroa, G.E. Romero and 
D.F. Torres, Mod. Phys. Lett. A 16, 984 (2001). 

[18] M.S. Morris and K.S. Thorne, Am. J. Phys. 56, 395 (1988). 

[19] Eq. (16.11) does not appear explicitly in Ref. m but it appears in a later paper m quoted 
in Ref. 1]. 

[20] J.F. Woodward, Found. Phys. Lett. 10, 153 (1997). 

[21] See Ref. B3 for an explicit discussion of the case m{t)a(t) =constant. 


14 


